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Abstract 

Additive genetic variances were estimated for several life history traits of north 
temperate Michigan and subtropical Florida populations of the freshwater copepod 
Mesocyclops edax (Crustacea: Copepoda). Estimates were obtained by rearing animals 
at 25°C and using a paternal half sib design. Considerable additive genetic variance 
was revealed in female and male body sizes of the Florida population and in female 
and male maturation times in the Michigan population. Female maturation time and 
clutch size in the Florida population and female body size and clutch size in the 
Michigan population exhibited low additive genetic variance. These findings provide 
further evidence of substantial additive genetic variance in life history traits of natural 
populations. The contrasting trends in additive genetic variance underlying life history 
traits of these populations are interpreted as reflecting differences in the potential to 
respond to selection and possibly also the different selective regimes of the Michigan 
and Florida environments. 


Introduction 

The variation in many life history traits of interest to evolutionary biologists is 
polygenic (Istock, 1981: Dingle and Hegmann, 1982; Lande, 1982) and suitable for 
quantitative genetic analyses. Plant and animal breeders have used such analyses to 
partition phenotypic variation into genetic and environmental components, and thus 
to predict responses to selection. Much experimental work on domesticated animals 
and plants, and on laboratory Drosophila populations, has documented changes in 
the proportion of additive genetic variance relative to other sources of variation after 
selection has been applied (see Falconer, 1980, for review). Often these findings bear 
out the predictions of Fishers’ Fundamental Theorem of Natural Selection (Fisher, 
1958): that traits exhibiting low additive genetic variance and thus showing low her- 
itability (ratio of additive genetic variance to total phenotypic variance) are presumed 
to have had much of their additive genetic variance trimmed away by directional 
selection. It is not known, however, whether this principle can be generally applied 
to natural populations. 

Some recent studies have uncovered substantial heritable variation in life history 
traits in natural populations and have attributed this variation to frequent changes in 
direction of selection (Dingle el al., 1977; Istock, 1981) or to weak selection and 
negative pleiotropy (Berven and Gill, 1983). As a population’s potential response to 
selection is. in theory, proportional to the heritable variation in fitness, the additive 
genetic variances of fitness characters are of interest when investigating life history 
( > evolution. If adaptation to a new environment results in changes in the amounts of 

additive genetic variance of particular traits, then changes in additive genetic variance 
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might be useful for inferring shifts in the relative contributions of traits to fitness. This 
line of reasoning was first applied to heritability analyses of a marine copepod by 
McLaren (1976) and of an estuarine copepod by Bradley (1978). Few organisms evolve 
rapidly in nature to provide an opportunity to observe such changes in additive genetic 
variance. An alternative is to examine geographically isolated populations of the same 
species that have presumably undergone evolution and genetic differentiation in re¬ 
sponse to their separate environments. This comparative approach attempts to provide 
insight into why specific traits might exhibit high additive genetic variance in one 
population and low additive genetic variance in another. 

In this study I estimate the additive genetic variances underlying phenotypic vari¬ 
ation in several life history traits in two populations of Mesocyclops edax( S. A. Forbes) 
(CrustaceaiCopepoda) from different locales. The two environments differ in several 
obvious ways: trophic status, thermal regime, and predator abundance and diversity. 
These populations show genetically based differences in the means of several traits 
(Wyngaard, 1986). Individuals from Florida exhibit more rapid maturation rates, larger 
body sizes, and larger clutch sizes than do those from Michigan. These traits are 
known to influence either survival and fecundity in natural zooplankton populations 
(Allan, 1976; Lynch, 1980; Kerfoot, 1980) and thus may be viewed as important 
fitness components. These two populations show no reproductive barriers in the lab¬ 
oratory, have fertile hybrids (Allan, 1984), and thus seem to be conspecific. The genetic 
differentiation in life history traits between the populations leads one to ask whether 
there have also been changes in the heritable variation of these traits. The purpose of 
this study was to estimate and compare the additive genetic variances of several fitness 
components among two geographically distinct copepod populations and to explore 
the usefulness of this information in making inferences regarding the potential to 
respond to selection as well as the past selective regimes experienced by these popu¬ 
lations. 


Materials and Methods 


Field study sites 

Individuals of Mesocyclops edax were collected from Lake Thonotosassa, Hills¬ 
borough County, Florida, on 10 July 1981 and from Douglas Lake, Cheyboygan 
County, Michigan, on 23-26 June 1981. Detailed descriptions of Lake Thonotosassa 
are given by Cowell et ai (1975) and Wyngaard et al. (1982) and descriptions of 
Douglas Lake are provided by Welch (1927), Tucker (1957), and Lind (1978). Here 
the Lake Thonotosassa population is designated as FT and the Douglas Lake population 
as MD. Plankton were shipped by air to the laboratory within 24 h of collection. 

Experimental design 

Additive genetic variance components and narrow sense heritabilities were esti¬ 
mated with a paternal half sib mating design (e.g., Falconer, 1980), using laboratory- 
reared animals that had undergone one generation of laboratory acclimation. Hap¬ 
hazardly sampled adult males and juvenile females from each population in nature 
were isolated in pairs and mated to yield the acclimation generation (G^. Sixty-nine 
MD (Douglas Lake, Michigan) and sixty-six FT (Lake Thonotosassa, Florida) matings 
were successful; four to twenty progeny from each of these matings were reared to 
adulthood. The adults of this acclimation generation were mated using a paternal half 
sib mating design in which two females (dams) were mated to a single male (sire). No 
adult was related to any other adult in the mating scheme. A total of 14 MD and 15 
FT half sib families were obtained. Twenty offspring from each dam were reared to 
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maturity and mated, taking care to mate only unrelated individuals. These offspring 
nsed (he second generation (G 2 ) on which the additive genetic variances of the 
ry traits were estimated. 

acclimation effects in maturation rate, body size (cephalothorax length), and 
h size were revealed between the two generations reared in the laboratory as mean 
values of these traits did not differ significantly (Wyngaard, 1983). Survival from 
hatching to adulthood of Gj individuals in both populations was 97%. Survival of G 2 
individuals of FT and MD populations was 96% and 91%, respectively. 

Laboratory rearing 

Upon hatching, progeny were replaced individually in petri dishes 20 mm deep 
and 60 mm in diameter containing artificial lake water (Wyngaard and Chinnappa, 
1982) and 10 4 cells/ml of Cryptomonas ozolini Skiijja (Starr Collection UTEX LB2194) 
isolated from Douglas Lake. The adult diet of M. edax was supplemented with the 
calanoid Diaptomus dorsalis (a common prey of M. edax in Lake Thonotosassa). 
D. dorsalis was mass cultured on Chylamydomonas reinhardti (Starr Collection No. 
89) and artificial lake water. Algae and crustaceans were grown at 25 ± 0.5°C and a 
light:dark cycle of 14:10. Individual M. edax were transferred to clean dishes with 
fresh medium and food every other day, and individuals were observed daily for 
survival and development. Upon maturity, unrelated males were placed with females 
until the first clutch was produced. Animals were preserved in 4% formaldehyde for 
body size and clutch size measurements. Cephalothorax length was measured to the 
nearest 0.003 mm at 40X with an ocular micrometer. 

Calculations of additive genetic variance components 

The sire variance component estimates the phenotypic covariance of half sibs and 
thus estimates one fourth of the additive genetic variance, because half sibs have one 
fourth of their genes in common, on the average (Falconer, 1980). Estimates from 
half sibs provide the best unbiased estimate of additive genetic variance. Dam com¬ 
ponent estimates, which include variance due to dominance (non-additive), maternal 
effects, and common environment, in addition to the phenotypic covariance of half 
sibs, were calculated for comparison with the half sib estimates. Because individuals 
were reared individually after hatching, common family environment is not likely to 
be important. Upon reaching maturity, males and females were paired; the males 
were removed when the females produced eggs. Cannibalism of males by females after 
mating occurred occasionally in the FI' matings and frequently in the MD matings, 
resulting in reduced numbers of males available for body size measurement. Sample 
size was too small to estimate the additive genetic variance of MD male body size. 
To minimize predation, some females and males were not mated but, instead, were 
preserved upon reaching maturity. 

As is commonly the case in such studies, class sizes were unequal; this resulted 
from variation in sex ratio and (to a small degree) differential survival of offspring. 
All data were used, and variance components for additive genetic variances were 
computed using maximum likelihood estimation (MLE) procedure VARCOMP of 
the computer package SAS (Hemmerle and Hartley, 1973; Helwig and Council, 1979). 
This method of extracting parameter estimates differs from that used in ANOVA but 
assumes the same underlying genetic model. MLE is preferable to ANOVA, because 
unlike ANOVA, MLE does not result in negative variance components, which are 
unrealistic. The superiority of MLE to ANOVA for use on unbalanced data is discussed 
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by Swallow and Monahan (1984). Standard errors of additive genetic variance com¬ 
ponents of half sibs were estimated according to Searle (1970). The significance levels 
of sire effects (or additive genetic variances) were tested by comparing the objective 
functions of the MLE analyses when the sire effects were estimated separately from 
the dam and within family variances to when they were combined with the dam 
variances. These objective functions closely approximate the mean squares (Hemmerele 
and Hartley, 1973), and their difference yields a Chi square value (x 2 ) which can be 
used to determine significance levels using 2 d.f. 1 present and compare estimates of 
additive genetic variance (V A ) in preference to heritability estimates (h 2 ), which are a 
ratio of the additive genetic variance (V A ) to the phenotypic variance (V P ). It is the 
genetic variance that is of interest in the present study; comparisons of V A are not 
complicated by the problems of comparing ratios. One trait, maturation time, was 
transformed to natural logarithms. 


Results 

The half sib (sire) components of variance on w hich the heritability analyses were 
based will be discussed first, as they provide the least biased estimates of additive 
genetic variance, and full sib (dam) components of variance will be discussed later. 
Neither population showed any significant heritable variation of clutch size or male 
maturation time (Tables I, III). There was a significant heritable variation for female 
maturation time in the MD population, but not in the FT population (Table I). In 


Table I 


Estimates of maximum likelihood variance components of maturation time (days) o/Mesocyclops edax 
from Lake Thonotosassa, Florida (FT) and Douglas Lake, Michigan (MD). 

Maturation time was transformed to natural logarithms 

Population 

Sex of 
offspring 

Source of 
variation 

df 

Variance 

component 

(SE) 

Percent 
of total 
variation 

x 2 

FT 

Female 

Among sires 

13 

0.0 

(0.001) 

0.0 

0 



Between dams. 

14 

0.0050 


39.1 




within sires 








Within families 

204 

0.0079 


60.9 




Total 

231 

0.0129 


100.0 


MD 

Female 

Among sires 

13 

0.0163 

(0.0077) 

42.8 

17.46*** 



Between dams. 

14 

0.0055 


14.5 




within sires 








Within families 

220 

0.0162 


42.7 




Total 

247 

0.0381 


100.0 


Fr 

Male 

Among sires 

13 

0.0 

(0.002) 

0.0 

0.01 



Between dams. 

14 

0.0077 


44.6 




within sires 








Within families 

231 

0.0096 


55.4 




Total 

258 

0.0173 


100.0 


MD 

Male 

Among sires 

10 

0.0103 

(0.010) 

28.0 

2.81 



Between dams. 

11 

0.0178 


48.3 




within sires 








Within families 

117 

0.0087 


23.6 




Total 

138 

0.0368 


100.0 



*** P < 0.001. 
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contrast, there was significant heritable variation for female body size in the FT pop- 
n v ji not in the MD population (Table II). The additive genetic variances of 
size could not be compared, because sample size was too small to obtain 
iate for the MD population (Table II). 

e degree to which the variance component among full sibs exceeds that from 
: half sibs indicates, in theory, the contribution of maternal common environment 
and non-additive genetic dominance effects to the observed variation. Comparison of 
estimates based upon sire and dam components indicated considerable nonadditive 
genetic and/or maternal components for FT male and female maturation time MD 
male maturation rate (Table I) and clutch size (Table III). The only instance in which 
the sire and dam compounds approximated one another was for MD female body 
size. Three traits exhibited a dam variance considerably less than that of the sire 
variance: MD female maturation time and FT male and female body size. 

Discussion 

The Florida (FT) population exhibited considerable additive genetic variance in 
body size, but not maturation time while the Michigan (MD) population showed the 
opposite trend. Both populations exhibited low additive genetic variance for clutch 
size. The standard errors of these variance components were large. However, in all 
cases (except maturation time of MD males) statistical tests of objective functions 
showed heritable variation only for traits that exhibited high amounts of additive 
genetic variance. Maturation time, body size, and clutch size influence survival and 
fecundity in natural populations and can be used to estimate the additive genetic 
variance of fitness components. 

Before I discuss these genetic variance estimates, their limitations should be noted. 
Presumably these laboratory derived estimates of additive genetic variance reflect those 


Table II 

Estimates of maximum likelihood variance components of body size (nm) ^/Mesocyclops edax 
from Lake Tlwnotosassa, Florida (FT) and Douglas Lake, Michigan (MD) 








Percent 



Sex of 

Source of 


Variance 


of total 


Population 

offspring 

variation 

df 

component 

(SE) 

variation 

X 2 

FT 

Female 

Among sires 

14 

31495 

(4472) 

13.7 

9.09* 



Between dams. 

15 

10409 


4.5 




within sires 
Within families 

160 

187377 


81.7 




Total 

189 

229281 


100.0 


MD 

Female 

Among sires 

13 

12210 

(4472) 

4.8 

0.97 



Between dams. 

14 

10906 


4.8 




within sires 
Within families 

174 

231608 


90.9 




Total 

201 

254742 


100.0 


FT 

Male 

Among sires 

12 

5849 

(1348) 

9.6 

6.74* 



Between dams. 

13 

0 


0.0 




within sires 
Within families 

125 

54892 


90.4 




Total 

150 

60741 


100.0 



* P <.05 
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Table 111 

Estitnaies of maximum likelihood variance components of clutch size o/Mesocyclops edax 
from Lake Thonotosassa, Florida (FT) and Douglas Lake, Michigan (MD). 


Population 

Sex of 
offspring 

Source of 
variation 

d.f 

Variance 

component 

(SE) 

Percent 
of total 
variation 

x 2 

FT 

Female 

Among sires 

13 

1.54 

(3.75) 

4.1 

0 



Between dams. 

14 

10.00 


26.4 




within sires 








Within families 

134 

26.25 


69.4 




Total 

161 

37.78 


100.0 


MD 

Female 

Among sires 

10 

0.00 

(3.54) 

0.0 

0 



Between dams. 

11 

5.58 


24.3 




within sires 








Within families 

74 

17.44 


75.7 




Total 

95 

23.02 


100.0 



of these populations in nature. The Florida population matures and reproduces over 
a wide temperature range (12-31°C) in nature, whereas the greater proportion of 
maturation and reproduction in the Michigan population occurs between 18-26°C. 
The lack of genotype-environment interactions, which would alter the additive genetic 
variances, is impossible to verify without repeating these experiments over a wide 
range of environmental conditions. Genotype-environment interactions have been 
reported for some life history traits in the marine copepod Eurytemora herdmani 
(McLaren, 1976) as well as for Drosophila melanogaster (Giesel et al., 1980, but see 
Rose, 1984), but were absent for some traits in D. pseudoobscura (Tantawy, 1961). 
General remarks regarding the relevance of genotype-environment interaction to her- 
itability studies are difficult to make, except that they can alter the interpretation of 
results and that they merit more investigation. 

The observation of traits exhibiting high additive genetic variance in one population 
and low additive genetic variance in another is intriguing in itself. The mean values 
of these traits (maturation time, body size, clutch size) show genetically based differences 
between populations (Wyngaard, 1986). The two populations interbreed readily in 
the laboratory (Allan, 1984); thus these populations comprise one species and com¬ 
parisons of their additive genetic variances are appropriate. 

The trends in additive genetic variances reported in the present study may be 
meaningful in one or several regards. First, the amount of additive genetic variance 
of a trait may reflect the intensity of selection in the immediate past upon that trait. 
Fisher’s Fundamental Theorem, which states that the rate of increase of mean pop¬ 
ulation fitness at any time is proportional to the additive genetic variance in fitness 
(Fisher, 1958), provides a context in which the relative amounts of additive genetic 
variance estimated for traits in this study can be compared and discussed. Certainly 
the large standard errors around these estimates of additive genetic variance make 
comparisons between populations tenuous; however, the differences between popu¬ 
lations in statistically detectable levels of additive genetic variance suggest such com¬ 
parisons are worthwhile. The Michigan and Florida populations inhabit lakes that 
differ in temperature regime, food level, and predator abundance (Wyngaard et al., 
1982; Elmore, 1983; Wyngaard, 1986). The relative contributions of maturation time, 
body size, and clutch size to fitness in each population may be related to these differences 
between the lakes. Confounding selection pressures may have molded these life histories 
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and comributed to obscuring any relationships that might exist among environmental 
variables nd genetic variances. Nevertheless, a plausible and testable scenario can be 
cor 'ii d that might explain, in part, the observed trends. 

ic possible explanation is that chronically low food levels in the Michigan lake 
a premium on reaching a certain minimum body size, maximizing energy de- 

i is of growth and maturation. Small clutches may be the result of a tradeoff between 
clutch size and egg size, favoring a few large young and enhanced survival in a low 
food environment. In contrast, rapid maturation and large clutches may be favored 
in the Florida population, in which predation is intense and mortality falls heaviest 
upon juveniles (Charlesworth, 1980; Wyngaard et al, 1982). The low additive genetic 
variances in maturation time and clutch size are consistent with this line of reasoning. 
At least, they suggest different potentials to respond to selection in the near future of 
microevolutionary time. Alternatively, low levels of additive genetic variance can result 
from random drift and founder effects. This latter explanation seems less likely, as 
plankton population sizes are typically large, AL edax is presumed to be highly vagile, 
and only small amounts of migration are necessary to overcome the effects of 
random drift. 

High levels of additive genetic variance pose problems of interpretation because 
they can be attributed to a greater variety of causes: weak selection, no selection, 
changing selection, negative pleiotropy, to name a few. Others have compared herita- 
bilities of various traits between populations in search of relationships between additive 
genetic variance and selection intensity. Berven and Gill (1983) argued that differences 
in heritabilities of larval duration and size at metamorphosis among three woodfrog 
populations suggested that selection intensities were very different in these populations. 
Arnold (1981) found similar heritabilities of feeding behavior in two geographically 
isolated snake populations. 

The trends exhibited by maternal effects and/or non-additive genetic effects in this 
study are troublesome but not unusual. The lower dam relative to the sire components 
of MD female maturation time and FT female body size in the present study may 
reflect non-random mating among parents with respect to genotype. If sires were 
mated with dams of similar genotype for a quantitative trait, then the genetic variance 
contributed by the dam would be masked in a nested design such as that used in this 
study. The results for AL edax probably represent a conservative estimate of the extent 
to which dominance and maternal effects are present in these traits. 

Perhaps the most important issue is whether the level of additive genetic variance 
reflects the potential response of the population to selection. The moderate and high 
heritabilities of female and male, respectively, temperature tolerance in the estuarine 
copepod Eurytemora affinis provided an opportunity to determine whether such her¬ 
itable variation was available to be acted upon by selection. Ketzner and Bradley 
(1982) exposed E. affinis to three different constant temperature regimes and to three 
variable temperature regimes over a two year period (20-70 generations) and dem¬ 
onstrated genetic divergence among the regimes as well as a decline in the additive 
genetic variance in one of the selected lines (Bradley and Ketzner, 1982). Other studies 
using selection experiments to address this question include 1 stock et al. (1975), Dingle 
et al. (1977), Doyle (1978), and Doyle and Hunte (1981). This growing body of em¬ 
pirical studies (review in Dingle and Hegmann, 1982) reveals considerable additive 
genetic variance to be present in natural populations of many species. The present 
study is consistent with these reports in that some of the life history traits examined 
here exhibited considerable additive genetic variance. The suggested contrasting pat¬ 
terns of additive genetic variance present in particular traits between the two AL edax 
populations indicate the need for additional studies of heritable variation of life history 
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traits in natural populations. Explanations of additive genetic variance are complex 
and will likely require both examination of genetic correlations and very carefully 
designed selection experiments as discussed by Istock (1982), Lande (1982), Rose 
(1982) and Lande and Arnold (1983). 
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